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ABSTRACT 

The last few months have seen the measurements of the radial velocities of all 
of the dwarf spheroidal companions to the Andromeda galaxy (M31) using the 
spectrographs (HIRES and LRIS) on the Keck Telescope. This paper summarises 
the data on the radial velocities and distances for all the companion galaxies 
and presents new dynamical modelling to estimate the mass of extended halo 
of M31. The best fit values for the total mass of M31 are ~ 7 - 10 x 10^^ M©, 
depending on the details of the modelling. The mass estimate is accompanied 
by considerable uncertainty caused by the smallncss of the dataset; for example, 
the upper bound on the total mass is ~ 24 x 10^^ Mq, while the lower bound 
is ~ 3 X 10^^ Mq. These values are less than the most recent estimates of the 
most likely mass of the Milky Way halo. Bearing in mind all the uncertainties, 
a fair conclusion is that the M31 halo is roughly as massive as that of the Milky 
Way halo. There is no dynamical evidence for the widely held behef that M31 
is more massive - it may even be less massive. 

Subject headings: galaxies: individual M31 - galaxies: kinematics and dynamics 
- Local Group - galaxies: structure - galaxies: halos 
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1. INTRODUCTION 

Dynamical studies of the two largest members of the Local Group - the Milky Way 
and the Andromeda (M31) galaxies - are of particular interest. The mass of M31 is well 
constrained on the scales of a few tens of kiloparsecs from the optical and H I rotation 
curves (Braun 1991) and from the kinematics of the planetary nebulae and globular cluster 
populations (e.g., Evans & Wilkinson 2000). To probe much larger radii and estimate the 
total mass of the M31 halo, it is to the satellite galaxy population that we must turn. Here, 
there have been dramatic developments in the last few years. The advent of wide-field 
imaging surveys (Armandroff, Davies, & Jacoby 1998; Karachentsev & Karachentseva 
1999; Armandroff, Jacoby, & Davies 1999) has led to the discovery of three new dwarf 
spheroidals (dSphs) - AndV, And VI (PegdSph) and CasdSph. In all, the entourage of 
galaxies accompanying M31 now numbers as many as fifteen, although the status of two of 
these (Pegasus and IC 1613) is perhaps unclear as opinions differ as to whether they are 
close enough to be bound (Mateo 1998; Grebel 1999; Courteau & van den Bergh 1999). 

For the Milky Way, the data set of the eleven companion galaxies, as well as the 
distant globular clusters, has often been exploited to estimate the total mass of the Milky 
Way halo. The older techniques of virial estimators have nowadays been superseded by 
more sophisticated Bayesian hkelihood estimators (Little & Tremaine 1987). The most 
recent analysis by Wilkinson & Evans (1999) found a total mass for the Milky Way halo 
of ~ 1911^ X 10^^ M©, with the large error bars reflecting the small sample size and the 
uncertainties in the proper motions. Hitherto, the satellites of M31 have received scant 
attention. Evans & Wilkinson (2000) looked at all the available published data on ten of 
the satellite galaxies and found a rather low mass of ~ 12.31^^ X 10" Mq. This raises the 
possibility that the Milky Way Galaxy may be the largest member of the Local Group. 

Here, we exploit recent stellar velocity measurements with the spectrographs on the 
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Keck Telescopes. The new observations are reported in detail by Guhathakurta et al. 
(2000a,b). The results enhance the dataset on the radial velocities of the companion 
galaxies of M31 by fifty per cent and are used to provide the most reliable estimate of the 
mass of the extended halo of M31 to date. 

2. DATA 

2.1. Radial Velocity Measurements 

The radial velocities of eight of the companions to M31 have been known for some 
time based on one or more of the following techniques: Hi 21cm measurements, optical 
absorption line spectra of the integrated light of galaxy cores/nuclei, bulges, or globular 
clusters, or optical emission line spectra of H ll regions. The data are recorded in the upper 
panel of Table |l|. These techniques are not easily applicable to galaxies of low surface 
brightness and/or low gas content. Consequently, until very recently, the radial velocities 
of the six dSph companions (And I, And II, And III, AndV, And VI, Cas dSph) remained 
unmeasured. They require 8-10 meter class telescopes and efficient spectrographs for 
spectroscopy of individual red giant stars. Cote et al. (1999) have recently provided the first 
kinematic study of And II. They used the High- Resolution Echelle Spectrometer (HIRES; 
Vogt et al. 1994) on the 10-meter Keck I Telescope to measure the radial velocities of a few 
red giants in each dwarf galaxy. Guhathakurta et al. (2000a,b) provide the data on the 
outstanding five dSphs. Using the Low- Resolution Imaging Spectrometer (LRIS; Oke et al. 
1995) on the Keck II Telescope, they obtained velocity measurements for nearly a hundred 
red giants in And I, And HI, AndV, and And VI. Guhathakurta et al. (2000a,b) have also 
carried out HIRES multislit echelle observations of 21 red giant branch stars in CasdSph. 
The radial velocities of the dSphs are recorded in the lower panel of Table |1]. 
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2.2. Distance Estimates 

In general, distance estimates for the M31 companions are derived from variable stars 
such as Cepheids and RR Lyrae, from the tip of the red giant branch (TRGB) method (Lee, 
Freedman, & Madore 1993), or from the brightness of horizontal branch (HB) stars (Lee, 
Demarque, & Zinn 1994). One has to rely on TRGB- or HB-based distances for the dSphs, 
since they lack the young stellar populations with which Cepheid variables are associated. 

Deep V and B exposures with Hubble Space Telescope/MVFPC2 have been used to 
clearly detect the HB in And I and And II (Da Costa et al. 1996, 2000). The distance 
estimates for these galaxies are the most accurate of all the M31 dSphs. The HB method 
provides greater distance accuracy than the TRGB method, but HB stars are faint and 
difficult to photometer accurately, especially in crowded fields, and this makes the method 
impractical for most ground-based telescopes. Moreover, the metallicity dependence of the 
HB absolute magnitude, My(HB), has been the subject of some debate (cf. Ajhar et al. 
1996). Distances based on the TRGB method are available for the rest of the M31 dSph 
satellites. Photometry of bright red giants has been carried out with the KPNO 4-meter 
telescope, 3.5-meter WlYN telescope, 6-meter BTA telescope at the Special Astrophysical 
Observatory, and with LRIS on Keck II for the following dSphs: And HI (Armandroff et al. 
1993), AndV (Armandroff, Davies, & Jacoby 1998; Grebel & Guhathakurta 1998), And VI 
(Grebel & Guhathakurta 1999; Tikhonov & Karachentsev 1999; Armandroff, Davies, & 
Jacoby 1999), and CasdSph (Grebel & Guhathakurta 1999; Tikhonov & Karachentsev 
1999). There are several hmitations to the TRGB method: it takes a large sample size 
to define the TRGB break accurately; since the method involves the location of a step in 
the stellar luminosity function, biases in the distance estimate are caused by blends and 
photometric errors; errors in statistical field subtraction can be a problem; super TRGB 
populations such as intermediate- age AGB stars can cause systematic errors. 
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For all the dSplis, there is substantial error on the distances (see Table |T]). Additional 
uncertainty is introduced through the poorly known distance of M31 itself, for which a 
variety of estimates is available using different methods; we adopt a distance of 770 kpc. 



3. MASS ESTIMATES 

Let us assume that the halo of M31 is spherically symmetric with the potential ip and 
density p given by 

The satellite galaxies are assumed to follow the same density law ps but with the scalelength 
Og = 250 kpc (cf., Evans & Wilkinson 2000). The velocity normalisation vq is chosen 
to reproduce the circular speed at 30 kpc of ~ 235 km s^^ (Braun 1991). We adopt 
two possibilities for the distribution of velocities. These must depend on the integrals of 
motion, namely the binding energy e and angular momentum I per unit mass. The first is a 
distribution function of form 

F{e,l) = l-'^f{e), (2) 

so that the orbital anisotropy (3 = 1— (vg)/ (v^.) is constant. This distribution ranges from 
radial anisotropy (/3 > 0) to tangential anisotropy {(3 < 0). The second is a distribution of 
Osipkov-Merritt form (see Binney & Tremaine 1987) 

F{eJ) = f{Q), Q = e-f/2Tl (3) 

where r^. is an anisotropy radius. This distribution is isotropic in the inner parts (r « r^) 
and radially anisotropic in the outer parts (r ^ r^). The functions /(e) and f{Q) are 
written out elsewhere (Wilkinson & Evans 1999; Evans & Wilkinson 2000). 

Mass estimation using Bayesian likelihood methods was introduced by Little & 
Tremaine (1987). Bayes' theorem states that the probability of any set of model parameters 
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given the data P(model| data) is 

P(model| data)P(data) = P(model)P(data| model). (4) 

Here, P(model) describes our prior beliefs as to the likelihood of the model parameters, 
while P(data| model) is the probability of the data given the model. For the 15 satellite 
galaxies, the data consist of three-dimensional positions rj with respect to the center of 
Andromeda, together with their heliocentric line of sight velocities VQ^i. These are corrected 
for the motion of the Sun around the Galactic Center and the infall of the Milky Way 
towards Andromeda. We assume a circular speed of 220 km s~^ at the Galactocentric 
radius of the sun (Pq = 8.0 kpc) and a solar peculiar velocity in km s"^ of (f/, V, W) 
= (—9, 12, 7). The line of sight velocity of M31, corrected for the motion of the sun, is 
^^r,M3i = —123 km s~^ (Courteau & van den Bergh 1999). This gives the velocities v^Q^i in 
the rest frame of the center of M31. Along any line of sight, let (ft, rj) be polar coordinates 
in the plane of the sky. Then the required probability is given by 

Pfdatalmodel) = — - — - / ' dvtVt / dr] f(ej) (5) 

The prior probabilities contain the a priori information about the model parameters. 
For the total mass M, we use P{M) oc 1/M^, while for the velocity anisotropy, we use 
P{/3) oc 1/(3 — 2(3)^ or P(ra) oc 1/ra (see Wilkinson & Evans (1999) for a discussion). 

The distances to the satellite galaxies are the main uncertainty. Our error convolution 
aims to take account of two factors. First, there is the quoted error estimate Asj associated 
with each distance given in Table |I|. Second, there is a small, but distinct, probability 
e that some of the published distance estimates are seriously in error due to systematic 
uncertainties; we assume that the probability of a rogue distance e is 0.1. Therefore, we 
choose our convolution function E to have the form 

E{z; i) = (1 - e)B{z; Si - Asi, Si + As,,) + eB{z; Smin, Smax)- (6) 
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Here, B is the hat-box function 




1 




'max 



'mill 



(7) 







otherwise, 



while Si is the pubhshed distance estimate of the ith satelhte with error Asj and Smin and 



Figure shows the resuhs for the constant anisotropy distributions. The most 
hkely mass is 7.0 x 10^^ Mq and the most hkely value of j3 is —0.95 corresponding to 
tangential anisotropy. An earlier analysis found that there is a tendency to underestimate 
the mass for small datasets (Evans & Wilkinson 2000). This typically leads to an 
underestimate of the order of a factor of two for a dataset of ten to twenty objects, and 
is combined with a spread that is of the order of fifty per cent. This gives a final result 
of ~ 7.0lll%^ X 10^^ Mq. The anisotropy is poorly constrained as the likelihood contours 
are distended horizontally. However, radial anisotropy is not favored as the Icr contour is 
concentrated in the tangentially anisotropic portion of the figure. Figure |^ shows the results 
for the Osipkov-Merritt distribution function. The most likely mass is ~ 9.7li\^g^ x 10^^ Mq 
and the most likely value of ra is ~ 680 kpc, suggesting the velocity distribution is roughly 
isotropic. Table ^ shows the effects of varying some of our modelling assumptions on the 
results for the mass. These include (i) the omission of Pegasus and IC 1613, which may 
not be bound to M31 (Courteau & van den Bergh 1999), (ii) the alteration of the satellite 
distances to those compiled by Mateo (1998), (iii) the increase of the error bars on the 
distances to 25%, (iv) the alteration of the scalelength of the satellite galaxy number 
density distribution and (v) the change of the velocity normalisation of the halo vq. The 
main change is the removal of Pegasus and IC 1613 from the dataset, which cause a roughly 
30% drop in the mass estimate. As can be seen from Table ^ the other alterations do not 
have a significant effect on the answer. 



-Smax are the maximum and minimum of the distance estimates listed in Table |I]. 
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4. CONCLUSIONS 

We have presented the first detailed study of the mass of the halo of the Andromeda 
galaxy that uses the radial velocity and distance measurements for all six dSph satellites 
(And I, And II, And III, AndV, And VI and CasdSph). This takes to fifteen the number 
of tracers at large radii and provides a reasonable sample with which to estimate the total 
mass. Such calculations have not been possible before, because of the unavailability of 
radial velocity data for the low surface brightness companions of Andromeda. Spectroscopy 
of individual red giant stars in these systems has only recently become feasible thanks to 
8-10 meter class telescopes such as Keck and state-of-the-art spectrographs (Cote et al. 
1999; Guhathakurta et al. 2000a,b). We reckon that the mass of the Andromeda halo is 
~ 7.0^3°/ X 10^^ Mq. This comes from analysing all the companion galaxies with a family 
of velocity distributions of constant orbital anisotropy. The error bars include the statistical 
effects of the small dataset, as well as the uncertainties in the distance estimators. There is 
some additional uncertainty from changes in the modelling assumptions but this is at the 
30% level at most. For example, altering the velocity distribution gives us a slightly higher 
estimate of ~ 9.7l^1^^ x lO^^ Mq. 

The most recent estimate of the mass of the Milky Way is ~ 19ti^ ^ 10^^ -^o (Wilkinson 
& Evans 1999). All the mass estimates inevitably come with substantial uncertainty, as 
they are inferred from the small datasets of satellite galaxies that are the only known probes 
of the distant halo. Bearing in mind such caveats, a fair conclusion is that the Andromeda 
halo is very roughly as massive as that of the Milky Way. There seems little evidence from 
the satellite galaxy motions for the widely held belief that the Andromeda halo is more 
massive than that of the Milky Way (e.g., Peebles 1996). Indeed, such evidence as there is 
points in the opposite direction, namely that the Andromeda halo may actually be slightly 
less massive. 
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Fig. 1. — This figure shows hkehhood contours in the plane of total mass M (in units of 
10^-^ Mq) and anisotropy parameter p. The contours are at heights 0.32, 0.10, 0.045 and 0.01 
of the peak. The most likely mass is 7.0 x 10^^ Mq and the most likely value of j3 is —0.90 
corresponding to tangential anisotropy. 
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Fig. 2. — This figure shows hkehhood contours in the plane of total mass M (in units of 
10^-^ Mq) and anisotropy radius ra- The contours are at heights 0.32, 0.10, 0.045 and 0.01 of 
the peak. The most likely mass is 9.4 x 10^^ Mq and the most likely value of ra is 680 kpc. 
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Table 1: Data on M31 and its companion galaxies. The heliocentric distances s are taken 
from Grebel (2000). The observed line of sight radial velocities are from Courteau & van 
den Bergh (1999; upper panel) or the new measurements discussed in the text (lower panel). 
Listed are Galactic coordinates {i,b), the hehocentric distances s in kpc, the observed line of 
sight radial velocities Vq in km s~^, the distances from the center of M31 r in kpc, corrected 
hne of sight velocities v^q (adjusted for the solar motion within the Milky Way and the radial 
motion towards M31), and object type. 



Name 


e 


b 


s 


Vq 


r 




Type 


M31 


121.2 


-21.6 


770±40 


-301±1 






Sbl-II 


M32 


121.1 


-22.0 


770±40 


-205±3 


5 


+95 


E2 


NGC 205 


120.7 


-21.1 


830±35 


-244±3 


61 


+58 


dSph 


NGC 147 


119.8 


-14.3 


755±35 


-193±3 


100 


+118 


dSph/dE5 


NGC 185 


120.8 


-14.5 


620±25 


-202±7 


173 


+107 


dSph/dE3 


M33 


133.6 


-31.5 


850±40 


-180±1 


225 


+72 


ScIl-III 


IC 10 


119.0 


-3.3 


660±65 


-344±5 


253 


-29 


dirr 


LGS3 


126.8 


-40.9 


810±60 


-286±4 


276 


-38 


dIrr/dSph 


Pegasus 


94.8 


-43.5 


760±100 


-182±2 


409 


+86 


dIrr/dSph 


IC 1613 


129.7 


-60.6 


715±35 


-232±5 


504 


-58 


IrrV 


Andl 


121.7 


-24.9 


790±30 


-380±2 


49 


-85 


dSph 


And II 


128.9 


-29.2 


680±25 


-188±3 


158 


+82 


dSph 


And III 


119.3 


-26.2 


760±70 


-355±10 


67 


-58 


dSph 


AndV 


126.2 


-15.1 


810±45 


-403±4 


118 


-107 


dSph 


And VI 


106.0 


-36.3 


775±35 


-354±3 


265 


-65 


dSph 


Gas dSph 


128.5 


-38.8 


760±70 


-307±2 


244 


-57 


dSph 



-17- 



Table 2: Illustration of the effects of changing some of the assumed model parameters. The 
parameter which has been changed is given in the first column. The second and third columns 
refer to distributions of the form /"^^/(e) while the fourth and fifth columns are for the f{Q) 
models. Also recorded are the most likely values of M (in units of 10^^ Mq), as well as the 
anisotropy parameters (3 and (in kpc). 



Comment 


Most likely 


Most hkely 
M 


Most likely 
re. 


Most likely 
M 


Canonical 


-0.90 


7.0 


680 


9.4 


Pegasus & IC 1613 


-0.95 


4.8 


224 


6.8 


omitted 










Sat. distances 


-0.85 


8.7 


659 


11.2 


altered 










As = 25% 


-0.90 


8.1 


610 


10.8 


Os = 150 kpc 


-0.95 


7.6 


627 


11.1 


= 270 km s"^ 


-0.95 


7.0 


697 


9.3 


= 200 km s-^ 


-0.95 


7.6 


648 


10.2 



